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We report on the first steps of an ongoing project to add gauge observables and gauge corrections 
to the well-studied strong coupling limit of staggered lattice QCD, which has been shown earlier 
to be amenable to numerical simulations by the worm algorithm in the chiral limit and at finite 
density. Here we show how to evaluate the expectation value of the Polyakov loop in the frame- 
work of the strong coupling limit at finite temperature, allowing to study confinement properties 
along with those of chiral symmetry breaking. We find the Polyakov loop to rise smoothly, thus 
signalling deconfinement. The non-analytic nature of the chiral phase transition is reflected in the 
derivative of the Polyakov loop. We also discuss how to construct an effective theory for non-zero 
lattice coupling, which is valid to <9(j3). 
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1. Introduction 

Due to the sign problem prohibiting lattice simulations at finite baryon density, the QCD phase 
diagram in the space of temperature T and chemical potential for baryon number /Xg remains largely 
unknown. Employing indirect methods like reweighting, Taylor expansions about jig = or sim- 
ulations at imaginary chemical potentials ji = i/J-i,^ € R, followed by analytic continuation, con- 
trolled calculations are only feasible as long as the quark chemical potential jj, = /j,b/3<T [1,2]. 
The latter two methods predict an at least initial weakening of the QCD quark hadron transition 
for quark chemical potentials ji < T [3-5]. A possible critical point at larger chemical potentials as 
well as a phase with net baryon number at low temperatures remain inaccessible to those methods. 
Furthermore, the chiral limit is inaccessible to direct simulations, thus obscuring its influence on 
the chiral properties of the QCD phase transition. These problems motivate the continued study of 
effective theories which are more amenable to finite density studies than full QCD. 

In order to control systematic errors, one would like to derive such effective theories directly 
from QCD, ideally with the possibility of systematic improvements by iteration of the approxima- 
tion scheme. Such a scheme is afforded by the strong coupling expansion. Recently, it has been 
used to derive a centre-symmetric 3d effective action for Yang-Mills theory, which provides a quan- 
titative description of the pure gauge deconfinement transition to within 10% of the corresponding 
full 4d simulations [6]. The extension to include Wilson fermions has also been reported at this 
conference [7], but is limited to heavy fermions. 

In this contribution we consider lattice QCD with staggered fermions in the strong coupling 
limit, i.e. for plaquette coupling /3 = 0. This is lattice QCD with one parameter chosen far away 
from its physical value. This is well motivated since the resulting model is confining, its massless 
limit exhibits chiral symmetry as well as its spontaneous breaking and, most importantly, it may be 
simulated at finite baryon densities by a suitable choice of algorithm. 

Starting from this model we can take into account gauge corrections, firstly with the aim 
of making contact with the successful description of the Yang-Mills action in terms of a strong 
coupling series [6], and secondly to gain insight into worm-inspired algorithmic methods [8] for 
Lattice QCD involving both, a fermionic and a gauge part in its action. Inclusion of fermionic 
effects and gauge corrections have also been pursued analytically in the mean field approximation 
[9]. Here, we aim to extend these investigations to a full evaluation of the effective theory. 

Starting point is lattice QCD with staggered quarks and the Wilson action, S = Sg + Sf, 

Sg = -f £ReTr£/ p (1.1) 

J p 

s f = l|l%W[zW^W^(^+A)-z(^+A)^WzW]+2^(x) z (x) I. (i.2) 

x { H J 

Expanding the partition function in a power series of j3 = 6/g 2 , 

Z = jDxDxDU(\-S G + 0(f3 2 ))e- s f, (1.3) 

the leading order expression corresponding to the strong coupling limit, /? = 0, is very simple due 
to the neglect of the gauge action. In this case there are only one-link integrals over the coupling 
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terms in the fermion action, which can be done exactly to produce a purely fermionic theory, 

Z(J8 =0) = j D X DxDUe Sf = f D X DxU I dU v (x)e- s " = f DzD%Uz(x,v) . (1.4) 



Its integrand now depends on colourless hadronic degrees of freedom only, mesons M x = XX( X ) 

3! 



and baryons B x = h&abcXaXbXc{x), 



3 

z(x,y) = £ Ok(M x ,M y y + aB x B y + pB y B x . (1.5) 
The final Grassmann integration can be done introducing a world line formulation to yield 

zo3=o)=£ n w*ii w *n w '> ^ 

{C}b=(x,v) x I 

with w x ,Wb,wi contributions from the chiral condensate, meson hopping and baryon world lines, 
respectively [10]. This so-called monomer-dimer model, representing the strong coupling limit of 
QCD, displays confinement as well as, in the massless limit for one flavour Nj ■ = 1, a U(l) chiral 
symmetry which spontaneously breaks in the vacuum. It is thus an interesting theory to investigate 
the finite temperature restoration of chiral symmetry. Since j6 = 0, the definition of a temperature 
scale requires introduction of an anisotropic lattice. The anisotropy parameter is defined in the 
weak coupling limit, y = a/a t . At strong coupling and in the mean field approximation one has 
instead y 2 w a/ a, and hence aT = y 2 /N t . The model in terms of baryonic degrees of freedom has a 
sign problem even at pL = 0, which however can be solved by an analytic reordering of terms [11]. 
For jU ^ 0, there is again a mild sign problem. Using a worm algorithm [8] allows for simulations of 
large lattices (spatial volume 16 3 ) at the relevant temperatures and densities [12]. Alternatively, the 



model Eq.(1.4) can be formulated and studied in continuous Euclidean time [13], taking the joint 
limit y,N t — > °°, aT = y 2 /N t fixed. There baryons become static and the sign problem vanishes. The 
strength of the sign problem at finite N t can be parametrised by quoting the average sign measured 
in the phase quenched theory, which can be expressed as a ratio of partition functions, 

Z V 2 

(sign) || = — = e ~T Af ^ 2 \ Z|| : phase quenched, A/~ 10~ 4 . (1.7) 

Using this approach, the phase diagram of the Nj = 1 theory in the chiral limit has been mapped 
out [12], revealing a second order chiral phase transition for small ji, which turns into a first order 
phase transition at a tricritical point, Fig.[j]. Note that, since there is the breaking/restoration of a 
symmetry involved, there is a true non-analytical phase transition separating the symmetric and 
broken phases. If one switches on small quark masses, chiral symmetry is broken explicitly and 
the second order phase transition changes into an analytical crossover. 

The long term goal of the present study is to extend the model Eq. ( |L~4] ) to finite values of the 
gauge couplings by including correction terms in the effective action, and possibly constrain the 
phase diagram in the weak coupling limit, as shown schematically in Figjl] (right). 
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Figure 1: Left: The phase diagram for Nf — 1 QCD in the strong coupling and chiral limit on an Nt = 4 
lattice [12]. Right: One possible extension for finite values of j3. For large j3 the phase diagram connects to 
that of the Nf — 4- theory in the continuum limit. 



2. Gauge observables in the strong coupling limit 

As a preparatory study towards including gauge corrections, we investigate the possibility of 
describing gauge field dynamics in the effective theory. In particular, one would like to study the 
confining behaviour of the theory across the chiral transition. However, in order to evaluate gauge 
observables in the strong coupling limit, we have to go back and write down their expectation 
values using the partition function in the strong coupling limit, but before gauge integration. Next, 
we rewrite the expression in terms of Zp , corresponding to the partition function in the strong 
coupling limit before fermion integration, 



DxDx / DUO[U]e 



Sf[U,x,x] 



JD X Dx {0)uZ F 



fDxD X Z F 

(O)u = — [ DUO[U]e- SF , Z F = [ DU e~ SF = \\z{x,v) . 



(2.1) 
(2.2) 



That is, for every gauge observable we wish to compute, the gauge integration has to be done with 
a different and now non-trivial integrand. Let us illustrate this with the Polyakov loop. We use a 
generating functional to define the expectation value by a functional derivative, 



(P) 
Zj 



d 

— InZj 
dJ J 



7=0 



Jd X Dx J D Ue- s ^ u ^ +Jpt+j1p ^ JdxDxZ f ((JP^ u + (J^P) l 



(2.3) 
(2.4) 



where we have approximated the exponential by its leading term, as we only consider infinites- 



imal /. After Grassmann-integration Zj takes on the form of Eq.( |1.6| ), this time with modified, 
/-dependent weights Wb, w x , Wi, and can again be evaluated using worm methods. Fig.|| shows the 
behaviour of (P) across the chiral phase transition as the temperature is varied. Interestingly, the 
Polyakov loop rises significantly, thus indicating deconfinement, even though we are in the strong 
coupling limit. However, the behaviour through the phase transition is smooth, similar to the find- 
ings by renormalisation group methods [14]. However, looking at the derivative of (P) in Fig.|| 
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Figure 2: Left: The Polyakov loop (P) as a function of temperature across the chiral transition for several 
volumes (am = 0). Right: Close-up of the transition region. 
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Figure 3: Left: The derivative d(P)/dT across the chiral transition. Right: The derivative of the internal 
energy de /dT shows a divergence consistent with a « 0. 



(left), we note that it develops a cusp. Thus, also the Polyakov loop signals non-analytic behaviour 
through the chiral phase transition, albeit in a higher derivative than the order parameter does. This 
can be understood by noting that the Polyakov loop can be written as a ratio of partition functions 
with and without a static quark Q, 

(P) ~ exp -{F Q -Fq)/T = ^-. (2.5) 

Since the free energy and its first derivative are continuous through a second order phase transition, 
so is the Polyakov loop. For comparison we also study the average number of meson hops in t- 
direction which is a measure for the internal energy e of the model. Its derivative Cy = de/dT, 
shown in Fig. || (right) as a function of temperature, is peaked at the transition with critical exponent 
a w as expected in the 3d 0(2)-universality class. 



3. Computing leading order gauge corrections 

We are now ready to address the calculation of gauge corrections to the strong coupling limit. 
To this end, we return to our initial expression, Eq. (1.3). For gauge observables in the strong 
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Figure 4: Left: Two ^(j3) -contributions to the effective action. Quark hops of the Dirac-part (arrows 
connecting circles) pair up with a plaquette-term to form colour-singlet degrees of freedom. Right: Average 
plaquette as a function of 1 / g 2 using the ^(/3)-contributionofthe{/(3) (purely mesonic) theory on a 8 3 x 16- 
lattice. For /3 — > they agree with the average plaquette obtained for SI/ (3) using HMC [15]. 



coupling limit it was sufficient to treat the relevant observables as source terms, c.f. the Polyakov 



loop, Eq. ([2.4]), expand the relevant term to first order in / and take the limit 7—^0 numerically. 
Here, we seek to take into account ^(j6) -corrections to the strong coupling action. We must hence 
study 



Z 



J DxDxZ F (zx V (-S G )) v « J DxDxZ Fe xp({-S G )u) , (3-D 



where in the approximation we used (exp(— Sg))u = ex P (L~=i e(( — ^G) k )u,c) and truncated at 
k = 1, i.e. to order 1 /g 2 . The implementation of gauge corrections now consists of two parts. First, 
we have to evaluate the gauge integrations present in the term ((— 5g)*){/,c which are non-trivial. 
Introducing the plaquette term leads to 19 new diagrams of the type shown in Fig.|] (left). After 
the gauge integration, they enter as new terms to the effective action. Note that this leading order 
correction has profound consequences on the physics of the effective theory. The gluon loop now 
allows for a separation of coloured quarks by one lattice spacing a. In particular, baryon hoppings 
on the lattice are now no longer self-avoiding, thus lending spatial structure to the baryon. 

The second part consists of the numerical evaluation of the partition function (or observables), 
now making use of the enlarged effective action. This is a non-trivial task, as it is not clear whether 
all additional terms can be regrouped in terms of the world line formulation that is then amenable 
to simulations by the worm algorithm. In our first attempt at a numerical evaluation, we have so 
far implemented only the mesonic contributions to the partition function, neglecting the baryons. 
In certain regimes of interest, in particular at zero net baryon density, this is a good approximation 
because baryons are heavy compared to mesons and therefore influence the dynamics only weakly. 

In Fig. |] we show our numerical results which, in the case of small but non- vanishing quark 
masses, we compare with those of full SU (3) HMC simulations, represented by the data points 
on the left. We see that for jS — > we obtain for the expectation value of the plaquette (U p )u — >• 
(U p )hmc> as expected. This is a quite non-trivial check, as these agreeing expectation values are 
simulated in terms of completely different degrees of freedom. The next steps to be taken are to im- 
plement the baryonic contributions, which will enable us to investigate finite coupling corrections 
to the phase diagram in the chiral limit. 
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4. Conclusions 

The strong coupling limit of QCD is interesting, because it represents a confining theory fea- 
turing a chiral phase transition in the massless limit, which is amenable to simulation by the worm 
algorithm, both at zero and finite baryon densities, even at low temperatures. We have shown that, 
even in the strong (i.e. infinite bare) gauge coupling limit, the chiral phase transition also implies a 
loss of confinement, as signalled by the rise of the Polyakov loop through the transition. Since the 
Polyakov loop is only indirectly associated with the order parameter of the chiral condensate, its 
expectation value goes through the transition smoothly, while its derivative develops a cusp. 

We are also on the way to including complete 0(j8)-corrections to the strong coupling limit, 
with a large number of analytic plaquette type contributions to the effective theory. Simulations 
using the worm algorithm including the mesonic part of the effective action give results consistent 
with full HMC simulations, and the implementation of the baryonic contributions is on the way. 
From an analytical point of view, higher order corrections are feasible. However, it remains to be 
seen whether the resulting additional terms to the effective action can be cast into a form amenable 
to be simulated by the worm algorithm, with a manageable sign problem. 

Acknowledgement: M.F. ,S.L. and O.P. are supported by the German BMBF, 06MS9150; W. U. is 
supported by the Swiss Nat. Fonds, 200020-122117. 
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